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Abstract 



Let the Lienard system u" + f(u)u' + g{u) = with an isolated 
periodic solution. This paper concerns the behavior of periodic solu- 
tions of Lienard system under small periodic perturbations. 
Key Words: perturbed systems, Lienard equation, polynomial sys- 
tems]]] 



1 Introduction 

Consider the second order differential equation of the type 

(E t ) x" + g(t,x,x',e) = 

where e > is a small parameter, g is a T-periodic function in t and 
g(t, x, 0, 0) = g(x) is independent of t. 

The existence problem of non constant periodic solutions of this equation in 
the case where g is independent on x' and is continuously differentiate has 
been studied by many authors. 

Indeed, in the latter case certain among them proved existence of solutions 
of x" + g(t, x, e) = 0. For a review see Chow-Hale [C-H] and Hale [H]. 
Example given by Hartman proved non existence cases of that equation if we 
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do not suppose g independent on x'. We then cannot generalize their result. 

Let the following equation, which is a perturbation of Lienard type 

(l e ) x" + f(x)x' + g(x) = eh( — ,x, x', e) 

where h is T-periodic in t, f and g are functions only dependent on x, sat- 
isfying conditions defined below. We look for periodic solutions of (l e ) for 
e small enough under some additional hypothesis. It is assumed that the 
unperturbed system has an isolated periodic solution. The perturbation is 
supposed to be controllably periodic in the Farkas sense [F], i.e. it is periodic 
with a period which can be chosen appropriately. 
We prove an existence theorem for this equation. 

Loud [L] already proves for f(x) = c, the existence of periodic solution 
of the equation 

x" + cx' + g(x) = eh(t), (1) 

where the perturbation does not depend on the state. He uses for that a 
variant of the implicit function theorem. More exactly, he considers a func- 
tion g(x) = xk(x) where k is continuously differentiable k(x) > 0, i^O. 
And either 

x— k(x) > 0, i^O 
ax 

or 

x— k(x) < 0, x^O 
ax 

always holds with the possible exception of isolated points. Let us notice that 
these conditions imply on one hand the monotonicity of the period function 
T for the system x" + g(x) — 0. If g is differentiable they imply on the other 
hand g"(0) = 0, that is a necessary condition. 
Let u(t) be a non-constant w-periodic solution of the equation 

x" + g(x) = 

and define _ 

F(s) = / u'(t + s)f(t)dt. 
Jo 

[L] observes that if for some so, F(so) = while F'(so) ^ then for 
sufficiently small e > there exists an cu-periodic solution v(t, e) of the 
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perturbed equation 

x" + g(x)=ef(t) = ef(t + u) 

2 Existence and non-existence of periodic solu- 
tions of E e 

2.1 A non existence result 

On the other hand, according to P. Hartman ([H], p. 39), equation (l e ) in 
general does not have a non constant periodic solution, even if xg(t, x, x') > 0. 
The following example given by Moser proves the non existence of a non 
constant periodic solution of 

x" + <j>(t, x : x') = 0. 

Let 

<t>(t,x,y) = x + x 3 + ef(t,x,y), e>0 

satisfying the following conditions for <fi G C 1 (i? 3 ), f(t+l,x,y) = f(t,x,y), 
with 

/(0,0,0) = 0, f(t,x,y)=0if xy = 



— — > oo when x — > oo 
x 



uniformly in (t,y) G R 2 , 



& f ^ f 

— > if xy > 0, and — = otherwise. 

by by 

x,y verifying | x \ < e, | y \< e. 

In fact, we have xf(t, x, y) and yf(t, x, y)) > if xy > 0, | x \< e, y arbitrary 
and = otherwise. 

Notice that ¥- is small. 

oy 

The function V = 2x 2 + x A + 2x' 2 satisfies V = — 4ex'f(t, x, x'), so that 
V < if xx' > 0, | x \< e and V = otherwise. 
Thus x cannot be periodic unless V — . 
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This example is significant because it shows in particular encountered 
difficulties in order to establish existence results of periodic solutions for ex- 
ample for the perturbed Lienard equation x" + f(x)x' + g(x) = eh(t, x, x', e). 
The period of the perturbed equation should be 'controlled' in order to state 
existence of periodic solution. 



2.2 Case where g is independent on x' 

Consider equation of the type 

x" + (j)(t, x, e) = 0, (2) 

where e > is a small parameter, is a continuous function, T— periodic 
on t such that <f>(t, x, 0) = g(x). 

More precisely, under the following hypotheses for the function g defined 

on R x (a, /3)x]0,e ]. 

' (1) is T — periodic on t 

< (2) <f>(t,x,G) = g(x) (3) 
k (3) if x ^ 0, we have g(x)x > 0. 

That means for e = the autonomous system 

{y'J-~g( X ) 

has the origin (0, 0) as a center. 
In using a version of fixed point theorem due to W. Ding, P. Buttazzoni and 
A. Fonda [B-F] proved there are periodic solutions of (2) provide that the 
period T of the autonomous associated system is monotone and e is small 
enough. is assumed to be (only) continuous 

More exactly, they show that if the function <fi(t, x, e) is continuous, then the 
periodic solutions of such equations may be located near to the solutions of 
the autonomous equation, provided that periodic solutions of (4) exist and 
the period function is strictly monotone. Moreover, there is a solution mak- 
ing exactly iV rotations around the origin in the time kT. 
Their results improve those of Loud [L] , that thought that the function 
had to be continuously differentiable. 

Moreover, under the light of the preceding example [H], it seems that the 
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methods described above does not generalize if one supposes dependent on 

x' : <f> = <j)(t, x, x', e). 

So, an other condition on the period appears to be necessary to obtain exis- 
tence of periodic solutions of the perturbed equation. 

One nevertheless can show an analogous result to the preceding one under 
more restrictive hypotheses. In particular on the appropriated choice of the 
period of a periodic solution of the perturbed Lienard equation. 



3 A controllably periodic perturbation 

One refers to a method due to Farkas inspired of the one of Poincare. The de- 
termination of controllably periodic perturbed solution. This method proved 
to be itself very effective particularly for the perturbations of various au- 
tonomous systems. Since we know it for example a good application for 
perturbed Van der Pol equations type [F-F]. We logically may expect that 
the Farkas method be again applied for perturbed Lienard equations. This 
has been considered and proved by Farkas himself [F] . 
Our proof we give here made it more simple and contains some modifications 
in using in particular methods of [Fl] to estimate existence regions of peri- 
odic solutions for that equation. 

The perturbation is supposed to be 'controllably periodic', i.e., it is periodic 
with a period which can be chosen appropriately. Under very mild conditions 
it is proved that to each small enough amplitude of the perturbation there 
belongs a one parameter family of periods such that the perturbed system 
has a unique periodic solution with this period. 

3.1 Basic hypotheses 

Let us consider the (unperturbed) Lienard equation 



In order to have an unique periodic solution we have to suppose the functions 
/ and g are of class C 2 . f is even, and g is odd. The integral 



(L) 



u" + f{u)u + g{u) = 0. 
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of / and g respectively are such that lim x ^ QO F{x) = oo, and lim x ^ OQ G(x) = 
oo. It is assumed that F has a unique zero. Then it is known [C-L] that 
(L) has a stable non constant periodic solution Uo(t) with period r . 
Equation (L) is usually studied by means of an equivalent plane system. The 
most used ones is : 



u' — v 



v' = -g(u) - f(u)v ^ 



and also 



u' = v- F(u) 

v> = -g(u) {b) 
In fact, there are equivalent to the 2-dimensional system 

{S) x = h(x) 
after introducing the notations x = col [xi,^] 

Xl = -u( t ) _ F ( u ( f )) 

x 2 = u(t) [n 

where x = col[xi,x 2 ] and h(x) = col[g(x 2 ), —X\ — F(x 2 (t)). 
Suppose 

u (0) = a, u'(0) = > 
so that the periodic solution of period To of the variational system 

y = h' x (p(t))y 

is 

p (t) = col[-u (t) - F(u (t)),u (t)\. 

We then have 

Pif) = col[g(u (t)),u (t)]. 
So, the initial conditions are 

p(0) = col[-F(a), a], p(0) = col[g(a), 0]. 
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Let us consider the following perturbed Lienard equation of the form 



where t ER, e ER is a small parameter, | e |< eo, r is a real parameter 
such that | r — r |< T\ for some < T\ < y. 
Moreover, the closed orbit 



belongs the region {(u, v) ER 2 : u 2 + v 2 < r 2 }. q is a function of class C 2 , 
r-periodic in t. 

By the same way as for (L), the 2-dimensional equivalent system to (Lr) is 



3.2 Existence of periodic solutions of (Lr) 

Now we will use Poincare method for the determination of the approximated 
solution of the perturbed equation (Lr). It is assumed the existence of the 
fondamental matrix solution of the first variational system of x = h(x) and 
the unique periodic solution p(t) corresponding to Uo(t). 
In order to get estimates for the existence of periodic solutions we have to 
calculate some constants. Following Farkas [F], the Jacobi matrix J has the 
following form 



/ = Id 2 and Y(t) is the fundamental solution matrix of the varational 
system with Y(0) = I 



(L R ) u + f(u)ii + g(u) = ey(-,u,ii) 



{(u,v) ER 2 : u(t) = u (t), v(t)=ii (t), tG[0,r ]} 



(S L ) x = h(x) + eq(-,x) 



where q = col\q u q 2 ], 




(8) 





(9) 
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It is proved, that if detJ(T ) ^ then there exist uniquely determined 
functions r(e, 0) and h(e, 0) defined in the neighborhood of (0, 0) such that 
the function 

u(f; (j),p + h(e,0),e,r(e,0)) 

is periodic solution of system (S£) so that r(0, 0) = To, h(0, 0) = 0. More- 
over, it is given an estimate for the region in which the the variables e and 
may vary . It is needed for that to evaluate the norm of the difference of 
Jacobi matrices J(e, 0, r, h) — J(0, 0, r , 0). 

A trivial r -periodic solution of (9) is col[g(u (t)), uo(t)]. A trivial calcu- 
lation gives the other linearly independent solution of (9) 

col[g(u (t))v(t),u0(t)v(t) + f^^m 

where 

v(t) = f[g Q {s)\- 2 g\u Q {t))exp[- f f(u Q (a))da]ds 
Jo Jo 

for t e [0,t ]. 

Then the fundamental solution matrix of (9) with Y(0) = I is 
«i g(a)g(u (t))v(t) 



Y(t) = 



9{a) 
9(a) 



^1 g(a)u (t)v(t) + g(a)^v(t) 



According to Liouville's formula the Wronskian determinant W(t) with 
W(o) = 1 is given by 

W(t) = exp[- f f(u (T))dr}. 
Jo 

The characteristic multipliers of (9) are p\ — \ and 

p2 = W(t ) = exp[- / f{u (r))dT]. 
Jo 

P2 < 1 if and only if 

f(u (r))dr > 0. (10) 
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The initial conditions give 
Thus, we get 



o P2 

j _ (9(a) 9 2 {a)v{r )\ 
\ P 2 ~ 1 J ' 
x (a) ^ 2 (a)w(r )(l -p 

o -(1-P2)- 1 



Therefore, 

|| J' 1 \\=2max |(?- 1 (a),(l-p 2 )- 1 ,( ? 2 (a)t;(r )(l-p 2 )- 1 ]. 
The inverse matrix of Y(t) is 

v . lfA _ ma ( 9(")*°(t)v(t) + 9(a)j^§v(t) -g(a)g(u (t))v(t) 
y \t) — W (t) I g(Mo(t)) 

Now we have to determine the constants for system (Lr). Following [F] 
let us denote 

S = {x= (x u x 2 ) eR 2 I x 2 2 + [-xx - F{x 2 {t))f < r 2 } 
g := max xeS | g(x 2 ) |, 

gi := max xeS I g'{x 2 ) |, (11) 
92 = max x&s \ g"(x 2 ) \ ■ 

fi := max x€S \ f(x 2 ) |, f 2 := max xeS \ f'(x 2 ) \ 



Qi 

L?2 



max xeS ,s<=R I g(s,a;) |, 

max xeS , s &i I I, (12) 



= max x€S ,seR I <?s( s ^) 



K := mai te[ ^ ]To] | F(t) |, := ma^^, | F x (t) 

Thus, we may deduce that 

P:=max te[Z z> To] | p(t) \< — , 
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The phase initial and the period r have to verify the following, which can 
be easely obtained from the above estimates, see [F] 

0<y, |t-t |<-. 

If in addition we suppose e and h are such that 

3 3 
-g | e | + ] h |< <7exp(--#iT ) 

(here a is the distance between the path of the periodic solution and the 
boundary of S) then a solution of (L R ) exists. 
We may resume in the following 

Theorem 1 If 1 is a simple characteristic multiplier of (9) (that 
means inequality (10) holds) then there are two functions r,h : U — > R 
and a constant t\ < y such that the solution u(t, 0, a + h(e, 0), e, r) o/ 

(L fl ) ii+ f(u)u + g{u) = ey(-,u,u) 

r 

exists for (e,<f>) G U, and satisfying properties \ r — r |< r 1; r(0,0) = 
r , //(().()) 0. 



3.3 Special cases 

As a corrolary we may deduce from the above some results about autonomous 
perturbations of the Lienard system 

(L Ra ) u + f(u)u + g(u) = 67 

where the perturbation is independent on the time variable 7 = j(u, ii, e, r). 
The equivalent plane system is of the form 

x = h(x) + eq( — , x) 

T 

where q = col[qi, #2]- 

Consider the following autonomous perturbed Lienard equation of the form 

(L Ra ) u + f(u)ii + g{u) = 67 = 7(u,u,e,r) 
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where t ER, e ER is a small parameter, | e |< e , r is a real parameter 
such that | r — t |< t\ for some < T\ < ^. 
Moreover, the closed orbit 

{(u,v) ER 2 : u(t) =u (t), v{t)=ii (t), tE[0,r ]} 

belongs the region S = {(u,v) ER 2 : u 2 + v 2 < r 2 }. 7 is a function of class 
C 2 . 

In this case the perturbation is independent on the initial phase 0. 
We then have the following 

Corrolary 2 Suppose inequality (10) to hold, then there are constants 
e and T\ < y such that to each e E [— e ,e ] there exist two functions 
r, h only dependent on e : r = r(e), h = h(e) such that the equation 

(L Ra ) u + f{u)u + g{u) = e7 = ey(u, it, e, r) 

has a unique periodic non constant solution u(t, e) with period r. Moreover, 
r(0) = t , and satisfy the properties \ r — To |< ri, t(0) = t , /i(0) = 0. 

On other hand, a second special case may occur when the perturbation 
does not depend on the state of the system. That means the perturbation is 
independent on u 

7 = 7(-)- 

T 

Then the above estimates can be easely calculated. 
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